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A linear time-varying aggregate traffic flow model can be used to develop Traffic Flow 
Management strategies based on optimization algorithms. However, there are no meth- 
ods available in the literature to translate these aggregate solutions into actions involving 
individual aircraft. This paper describes and implements a computationally efficient disag- 
gregation algorithm, which converts an aggregate (flow-based) solution to a flight-specific 
control action. Numerical results generated by the optimization method and the disaggre- 
gation algorithm are presented and illustrated by applying them to generate TFM schedules 
for a typical day in the U.S. National Airspace System. The results show that the disag- 
gregation algorithm generates control actions for individual flights while keeping the air 
traffic behavior very close to the optimal solution. 


N omenclat ur e 

number of centers 

directional link from center i to center j 
time horizon of simulation and optimization 

time interval when weather impacts traffic flow between center i and center j 
initial condition 

maximum number of aircraft allowed in center i at time k 
maximum number of departures allowed to enter center i at time k 
number of departures scheduled to enter center i at time k 
maximum number of arrivals allowed to leave center i at time k 
total number of scheduled arrivals in center i 
maximum flows allowed between centers i and j at time k 

minimum time that a flight has to spend in center i if the flight goes through it 
optimal number of aircraft in center i at time k 
equals f3u(k)xi(k), optimal number of arrivals in center i at time k 

I. Introduction 

Today, air traffic flow prediction is done by propagating the trajectories of the proposed flights forward 
in time and using them to count the number of aircraft in a region of the airspace. Examples of systems that 
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use this physics-based modeling approach for demand forecasting include the Center TRACON Automa- 
tion System (CTAS), 1 the Future ATM Concepts Evaluation Tool (FACET), 2 and the Collaborative Routing 
Coordination Tool (CRCT). 3 The accuracy of these predictions is impacted by departure and weather uncer- 
tainties. 4,5 These trajectory-based models predict the behavior of the National Airspace System adequately 
for short durations of up to 20 minutes. With the short prediction accuracy, it is difficult to make sound 
strategic decisions on air traffic management. 

Aggregate models simplify the analysis and design of many complex systems. The development of ag- 
gregate flow models for traffic flow management has been the subject of considerable interest, since the first 
model 6 appeared in the literature. These models are based on aggregation of traffic flows and simulate the 
behavior of volumes of flights. Compared with trajectory-based models, aggregate flow models are compu- 
tationally efficient, and their complexity does not increase with the number of aircraft. The aggregate flow 
models previously published 7 represent and predict the traffic behavior to a high degree of accuracy and can 
be tailored to the time-scales and regions of interest. A comparison of the characteristics of the different flow 
models can be found in previous studies. 8,9 The stability and response characteristics of the aggregate flow 
models are presented in the article by Chatterji and Sridhar. 10 The aggregation in flow models generally 
results in the loss of information about the route structure of individual aircraft. Additional research is 
needed to translate the aggregate control policies into actual traffic flow management flight planning deci- 
sions involving aircraft departure times and routes. Currently, there is limited experience in the application 
of aggregate and reduced order flow models for generating flight-level traffic flow management decisions. 

In this paper, an aggregate traffic flow model developed in earlier work 7 is implemented. The system 
parameters of the model are identified using historical air traffic data. Using filed flight plans as an input 
to the model, the traffic situation is predicted and compared with real data. A two-step approach is used 
to solve traffic flow management problems. The first step is solving an optimization problem, in which the 
system dynamics of the aggregate model are part of the constraints. In the second step, a computationally 
efficient disaggregation algorithm is designed to convert flow-based optimal solutions from the first step to 
flight-specific control actions. 

This paper is organized as follows. The next section briefly introduces the aggregate flow model. The 
optimization algorithm and the disaggregation algorithm are described in Section III and Section IV, respec- 
tively. Section V presents application results of the optimization framework. Finally, concluding remarks 
and future work are discussed in Section VI. 

II. An aggregate flow model 


A. A dynamic system model 

An aggregate traffic flow model was developed in an earlier work. 7 It is a linear-time-variant dynamic system 
model (LDSM). The number of aircraft at different times in each center is represented by a state variable.* 1 
The number of landings in a center and transitions from the center to the neighboring centers in an interval 
of time AT are assumed to be proportional to the number of aircraft in the center at the beginning of the 
interval. Using the principle of conservation of flow in a center, the number of aircraft in center i at the next 
instant of time k + 1 can be related to the number of aircraft in i at k via the difference in the number of 
aircraft that came into the center and the number of aircraft that left the center as follows: 

N N 

Xi(k+ 1 ) = Xi(k) -^2pij(k)xi(k) + ^ (3ji(k)xj(k) + di(k). ( 1 ) 

3 = 1 3 = 1 , 3 ^* 

where fiij{k) is the fraction of the aircraft in center i at the start of a time step that leave center i and 
enter center j during the k - th time interval AT. The number of arrivals in center i during the k-ih time 
interval, denoted by f3a(k) ■ Xi(k ), is a fraction of Xi(k). The departure within center i is denoted by di(k ), 
which is independent on Xi(k). For simplicity of illustration, “during the k- th time interval AT” will also 
be understood as “at time fc” in this paper. 

a The LDSM can be used to model different levels (sizes) of airspace, e.g., sector-level LDSM and user- defined- level LDSM, 
with respect to different levels of objectives and interests. The optimization algorithm desigined in this paper can be adpated 
for different LDSM models accordingly. 
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B. Implementation and validation of LDSM 

The inputs to the LDSM include departures di(k) within center i at time k and the fractions f3ij(k). The 
direct output from LDSM is the aircraft count in the centers at each time step. It is also straightforward to 
generate other outputs, such as inter-center traffic flows, number of arrivals, etc., based on the information 
of aji(fc), pij (k) and di(k). 

To implement the LDSM, the fractions faj ( k ) are obtained as the aggregated fractions of aircraft going 
from center i to j during the same k-ih time interval in each day, using historical air traffic data. In this 
study, the departures di(k) are computed from filed flight plans (deterministic). This is different from the 
original model in Ref. 7, which includes both deterministic and stochastic components. The LDSM model 
is implemented in a deterministic manner because a deterministic optimization method (Section III) will be 
used for Traffic Flow Management (TFM). The LDSM is implemented in C++; all the inputs were generated 
by processing flight data using FACET. 2 All 22 centers of the United States airspace, oceanic centers and 
part of the Canadian centers are included in the LDSM. 

The LDSM is validated against real air traffic data. Figure 1 shows a comparison between real and 
predicted aircraft counts, with two time intervals AT = 1 minute and AT = 15 minutes, in Oakland center 
(ZOA) from PDT 5:00PM on August 23, 2005, to PDT 4:59PM on August 24, 2005. For this validation, the 
fractions ffij (k) are obtained using the historical traffic data from September 6, 2001, and the departures 
di(k) are generated using the flight data from August 24, 2005. It is shown that the LDSM correctly predicts 
the trends of the evolution of aircraft counts in each center, and the average relative error between prediction 
and real data is less than 2%. 

It took less than 30 seconds to process one-day of historical data using FACET and to compute (3 with the 
C++ code on a 1.8GHz CPU, 2GB RAM IBM ThinkPad T42 laptop, operating system Linux. Performing 
a prediction of air traffic for a whole day (24 hours) took approximately five seconds. Compared with most 
trajectory-based models that predict the behavior of the National Airspace System (NAS) adequately for 
shorter durations, the LDSM is efficient to predict the traffic situation for a longer time horizon, which is 
helpful to make strategic decisions on air traffic management. 

The aggregate flow model can be used to design traffic flow management strategies to minimize different 
objectives subject to various constraints. Several optimization methods have been designed for traffic flow 
management, 11-15 most of which are based on networked air traffic flows. In the next section, a new 
optimization method will be proposed, which fully utilizes the aggregate flow properties of LDSM, and is 
suitable for future disaggregation to generate flight-specific control. 

III. Optimization method 

In this section, using the aggregated flow model LDSM, an optimization method will be designed to 
optimize the parameters /3 and d{(k) in equation (1) for traffic flow management. 

A TFM problem is formulated as follows. 

min J(x) (2) 


subject to 

Initial condition (Section III-A.l) 

LDSM dynamics (Section III-A.2) 

Center capacity constraint (Section III-A.3) 

Inter-center flow constraint (Section III-A.4) 

Departure constraints (Section III- A. 5) 

Arrival constraints (Section III-A.6) 

Minimum flight time constraints (Section III- A. 7) 

The constraints and the cost function (2) are described in the following subsections. 
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Figure 1. Oakland center traffic. 
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A. Constraints 


1. Initial condition 
The initial condition 

z(0) = x 0 (3) 

represents the number of aircraft in each center at the beginning ( k = 0) of the optimization. 

2. LDSM dynamics 

The LDSM is subject to the following dynamics: 

N N 

Xi(k + 1) =Xi(k) ~'Y^pij{k)x i {k) + ^ 0ji(k)xj(k) + di(k), he {(),••• ,T},ie {1, — ,N} (4) 

3 = 1 J = 

which is presented in Section II. 

3. Center capacity constraint 

The time-varying center capacity constraints are formulated as 

0 <®<(fc) <Cl(fc), fce{0,.-- ,N} (5) 

which restricts the number of aircraft in a center to be below the maximum number of aircraft allowed. 

4 . Inter-center flow constraint 
Inter-center flow constraints are formulated by 

PiA^Xiik) + 0ji(k)xj(k) < k € {7y, • • • , 2* } (6) 

which forces the flow between centers i and j below the maximum amount of flow allowed during the time 
interval • • ,2$]- 

5. Departure constraint 
Departure-related constraints are as follows: 

0<di(k) < Di(k), ke {0, - ,T},ie {1, - ,AT} 

t t 

^(^<^2 Si(k), i € {1, • • • , N}, t € {0, • • • , T - 1} 

k= 0 k= 0 

T T 

J2di(k) = J2 S i( k )’ 

fc=o k= 0 

where equation (7) restricts the number of departures below the departure capacity at each time step; 
equation (8) enforces that the cumulative number of departures in a center by time t cannot exceed the 
cumulative number of scheduled departures by that time; equation (9) enforces all flights to depart. 

6. Arrival constraint 
In the arrival constraints 

0 < Oi(fc) < Ai(fc), k e {0, • • • , T},i e {1, ■ ■ ■ , N} (10) 

T 

J2a i (k) = R i , i € {1, • • • , N}, (11) 

fc=0 

equation (10) restricts the number of arrivals below the arrival capacity, and equation (11) enforces total 
number of arrivals by time T (left-hand side of the equation) to be equal to the total number of scheduled 
arrivals (right-hand side of the equation). 
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7. Minimum flight time constraint 

This constraint requires a minimum time that a flight has to spend when flying through a center: 

N 


^ ^ fiij (k)xi (fc) — 0, aflk) — 0, 

3 = 1 
N 

^ v = V,di{k) = 0, 

3 = 1 

t+Ti / N \ 

E ^ 2 / 3 ij(k)xi(k) + ai(k) 

k=T t \j = 1 ) 


fee {(),••• ,Ti-i},ie {l,--- ,N} 


(12) 

k G {T — Ti + !,•••, T}, *€{!,••• 


(13) 

t / N \ 

^E E ’ 

fc=0 \j= 1 I 

t e {o, ■ ■ 

■ ,T-Ti},i = N} 


(14) 


Details about this constraint are presented in Appendix B. 

Remarks 

In all constraints, the fraction variables flij(k) always appear together with Xi(k). Therefore alternative 
variables f%j(k) are defined where fij(k) = (3ij(k)xi(k); fij(k) is actually the amount of traffic flow that 
leaves center i and enters center j at time k. Replacing /3ij(k)xi(k) by fij{k), all constraints (3)-(14) are 
linear in £*(&), di(k) and fij(k). 

B. Cost function 

The cost function J(x) in (2) can be defined in various forms according to specific objectives. For example, 
to minimize the difference between center counts of flights and a desired center count profile, it can be defined 

as nt 

^(*) = EE II** lla> ( 15 ) 

i = 1 k = 0 

where xf(k) is the desired center count profile. 

Another example of a cost function is 


N T 


J( x ) = E £**(*)• 

i= 1 k—0 


which is equivalent to minimizing the total flight time (see Appendix). 

When minimization of the departure delays is also part of the objective, the cost function can be formu- 
lated as 

j(*) =EE*iW+EE f E ^i(^) -E*( fc ) 

i = 1 k—0 i = 1 t=0 \k—0 k—0 


Jo(x) 


Ji(d) 


where Jq(x) is the total flight time, and Ji(d) is the total departure delay. 

Delays in different centers can also include designated weights to impose some requirements on the equity 
with which delays are distributed across centers. 16 The objective function can be formed as follows, 


^ 0*0 = EE< nroute (k)xi(k) + tt(t ^departure (5 . (fc) _ d . {k)) 

i= 1 k—0 i= 1 t—0 \k = 0 

where wf n route and wf eparture are the weights for en route delay and ground delay in center i, respectively. 

The objective functions above are subject to the same constraints described in the previous subsection. 
The first objective function in equation (15) is convex (quadratic), and the optimization problem is convex. 
The other two objective functions are linear, and the optimization problems are linear. 
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IV. D isaggregat ion algor it hm 


The LDSM and the optimization model are formulated based on aggregated traffic flows. An individual 
flight’s identity is not preserved. A disaggregation algorithm, which converts an optimal aggregate flow- 
based solution ( d and /3) b to a flight-specific control action, is presented in this section. The optimization 
algorithm outputs for each center the number of departures, the number of aircraft going to each neighbor, 
and the number of aircraft that should be in it, during each time period. This algorithm then determines 
which flights should depart from a center during each time period, the route (a sequence of centers) and the 
transition times for each flight. 

The disaggregation algorithm includes the following steps: 

1. Round di(k) and fij(k) = (3ijXi(k) to integers. 

2. Generate a path map: find m possible flight paths (routes) for each origin-destination center pair from 
historical air traffic data. Rank the paths from “shortest” to “longest” based on the average flight 
times. In this paper, m = 3. 

3. Generate a departure queue for each center by picking Y^k = o <&(&) flights taking off from center i. The 
flights with earlier scheduled departure time are put at the front of the departure queue. Flights on 
top of the departure queue have higher priority to take off. 

4. At each k for each i: 

(a) Take out di(k) flights from each departure queue and push these flights into the center flight queue 
(a queue of flights in a center) from which the flights depart. 

(b) Disaggregate inter-center flows (pseudo code in Algorithm 1): 

i. For each center flight queue, pop out the flights which land in this center; push these flights 
in the arrival flight queue of this center. 

ii. For each inter-center flow fij(k): search flights in the center flight queue of i, find possible 
flights that may go through center j. Based on the path map built in Step 2, rank these flights 
in an order such that flights with a shorter path to destination will have a higher priority to 
enter center j. 

iii. Pop out fij(k) flights from center i flight queue and push them into center j flight queue. 

(c) For each arrival queue, check if a flight is eligible to land (current criteria: it has to be in the 
arrival queue of center i for at least X- units of time). If eligible, pop out this flight from the 
arrival queue and let it land. 

5. A “clean-up” process to remedy the numerical rounding problem. 0 For each airborne flight fid in 
center i, 

(a) Compute the earliest time t l j id when fid can leave center i based on its entry time and the 
minimum flight time requirement. 

(b) If fid is eligible to land in center let it land. 

(c) Otherwise, find a minimum cost path from center i to fid's destination. The “cost” includes 
minimum flight times in the centers along the path, and the marginal capacities of the centers to 
accept additional flights, i.e., Cj(k) — Xj(k) for center j at time k, where Xj(k) is the number of 
disaggregated aircraft in Step 4. If Cj(k) — Xj(k) < 0, a very high cost is assigned for center j at 
time k. 


b This disaggregation algorithm works for solutions other than optimization as well: given the number of departures and 
inter-center transition flows, the algorithm converts the flows into flight- specific actions, if possible. 

c As a result of numerical rounding to di(k) and fij{k ), center capacity constraint of equation (5) may be violated, and some 
flights may be airborne at the end of optimization time horizon. 
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for each flight fid in center queue C enter Q(i) do 
if dest(fid) = i then 
CenterQ (i).pop(fid) 

ArrivalQ(i).push(fid) 

else 

if path(j, dest(fid)) exists then 
tmpSet(i,j) = tmpSet(i,j) U {fid} 

end if 
end if 
end for 

rank flights in tmpSet(i, j)\ 

if length of path(j, dest(fidi)) < length of path(j, desttfidz)) then 
priority (fidi) > priority (fid^) 

end if 

for n = 1 to fij(k) do 

fid n = highest priority flight in tmpSet(i, j) 

CenterQ(j ) .push(fid n ) 

C enter Q(i).pop(fid n ) 
tmpSet(i,j) = tmpSet(i, j)\{fid n } 
end for 

Algorithm 1: Disaggregation for inter-center flows. 


V. Application to TFM 

The optimization model in Section III and the disaggregation algorithm in Section IV are implemented 
on the same platform as in Section II. B. In the implementation of the optimization model, ILOG CPLEX 
Concert optimization library 17 is called in C++ to solve the optimization problem. Several application 
examples of the model and the algorithm are described in the rest of this section. Four-hour TFM problems 
for the airspace including 20 continental centers of the United States were solved in the examples. The 
states of the LDSM are updated every 15 minutes. All 8422 flights that departed from the continental 
centers during the time horizon were included in the study. 

A. Minimization of delays with center capacity constraints 

In the example presented in this section, the objective is to minimize the total flight time and the departure 
delays under constraints of center capacities. The cost function is defined as follows: 

j(x) = £X>(fc) +££ (l>(*) - £*(*)) ( 16 ) 

i= 1 k= 0 i= 1 t=0 \k= 0 k= 0 / 

Figures 2 and 3 show four center count profiles in Cleveland (ZOB) and New York (ZNY) centers: 

1. The profile labeled Real is the recorded real center aircraft counts. 

2. The profile labeled Simulated is a prediction of center aircraft counts using LDSM, i.e., the fractions (3 
in (1) are computed from historical data, and the departures di(k) = Si(fc), where Si(k) are scheduled 
departures. 

3. The Optimized profile is an output from the optimization with the objective function of equation (16). 
Notice that the optimized center counts are below 200, which is the user-defined center capacity of 
ZOB in the example. The simulated and real center counts exceed the center capacity during the first 
two units of time. 

4. The Disaggregated profile is the output from the disaggregation algorithm, which generates flight- 
specific action. 


8 of 16 


American Institute of Aeronautics and Astronautics 




ZOB 



Figure 2. Comparison of center count profiles in Cleveland center. 


ZNY 



Figure 3. Comparison of center count profiles in New York center. 


Table 1 summarizes the costs defined by equation (16) for the real, simulated, optimized and disaggregated 
data. Clearly, the optimized solution has the minimal cost among the four, while the disaggregated profile 
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Optimization run time 



Figure 4. Average computational time for optimization. Horizontal axis denotes 31 days in May 2005. 


Disaggregation run time 



Figure 5. Average computational time for disaggregation. Horizontal axis denotes 31 days in May 2005. 
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Center counts (z) 

Real 

Simulated 

Optimized 

Disaggregated 

J(x) 

1.28 x 10 7 

1.3 x 10 7 

9.1 x 10 6 

1.0 x 10 7 


Table 1. Comparison of the costs with different center count profiles. 


obtained with flight-specific information has a larger cost. The simulated data is very close to the real data 
because the LDSM predicts flows and aircraft counts in centers very close to the real situation. Depending 
on the values of TVs in the minimum flight time constraint in equation (14), the optimal solution has a lower 
cost compared with the real data and simulated data in general. 

The average computational time for optimizations of several test cases in each day of May 2005 is shown 
in Fig. 4; the average computational time for disaggregation is in Fig. 5. The computations are performed 
on a 1.8GHz CPU, 2GB RAM IBM ThinkPad T42 laptop with Linux operating system. 

B. Optimal flow routing 

Figure 6 displays a weather scenario in which Chicago (ZAU), Indianapolis (ZID), Kansas City (ZKC) and 
Memphis (ZME) centers are impacted by the weather, and their capacities are reduced correspondingly. 
Using LDSM, it is observed that there is a major flow going through ZKC-ZAU-ZOB under normal weather 
conditions. Due to the weather impact, the center capacity of ZAU is assumed to have a 75% reduction. 18 
Under this assumption, the optimization model and the disaggregation algorithm generate a routing strategy 
such that around 70% of the flights, which use the ZKC-ZAU-ZOB route under normal weather conditions, 
used a route ZKC-ZMP-ZOB instead. Figure 7 demonstrates this scenario: the blue dashed line represents 
the route under normal weather condition, while the red solid line represents the route under weather impact. 



Figure 6. An example of NCWD Weather Report. 


VI. Conclusion and future work 

In this paper, a linear time-varying aggregate traffic flow model is implemented and validated against 
real air traffic data. Using the aggregate flow model, an optimization model is proposed for traffic flow 
management. A disaggregation algorithm is developed to convert a flow-based solution to a flight-specific 
solution. Disaggregation is an important part for aggregate modeling for implementable control actions. The 
optimization model and the disaggregation algorithm are implemented and applied to TFM problems. 

Future work will focus on the following items: (^) The aggregate model in this paper only includes the 
deterministic components in the original model, and the corresponding optimization is also deterministic. 
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Figure 7. Dashed line (blue) is the original route. Solid line (red) is a reroute. 


Future work will include both deterministic and stochastic components in the aggregate model, and stochastic 
optimization algorithms will be explored for generating traffic flow management controls. {%%) The current 
disaggregation algorithm does not explicitly include capacity constraints (center capacity and flow capacity). 
A refined disaggregation algorithm will be designed in which the capacity constraints will be respected. 

Appendix 


A. State and travel time 

One of the characteristics of the LDSM model is that the states of the dynamical system are directly related 
to flight time. State X{(k) represents the number of aircraft in the center at time k. Assuming the flights do 
not transfer to a neighbor center or land, Xi(k) is also the number of aircraft in center i during a left-closed 
and right-open time interval [fcAT, (fc+ 1)AT). Therefore, Xi(k)AT is the flight time for all aircraft in center 
i in the k-ih time interval. Clearly, the total flight time spent in center i from time T-° to T- 1 is: 

n 

AT Y 

k=T° 

For all the centers, the total flight time (during the entire time horizon) is 

A rf)f>(k). 

i—1 k = 0 

Therefore, minimizing the following cost function 

j( x ) =YYj Xi ^ 

i= 1 k = 0 

is equivalent to minimizing the total flight time. 

B. Minimum flight time constraint 

Because every aircraft flies at a reasonable speed (not too fast) in a center, the optimization method enforces 
that each flight has to spend a minimum amount of time in each center it passes. 
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In the constraint (12), J2f=i Pij(k)xi(k) = 0, k G {0, • • • , T* — 1} means that the outflows from center i 
(the left-hand side of the equation) in the first T* units of time is zero, because any flight in center i has to 
stay in the center for at least Ti units of time and cannot exit from the center before time T-. Combining 
the non- negativity constraints on (3 and x , this constraint is identical to the following constraint: 

pij c k)xi(k ) =0, k e {0, . . ■ , Ti — 1}, j = {1, . . . , N}. 

The second part of (12), di(k) = 0, k G {0, ■ ■ ■ , Ti — 1} means that there are no arrivals in center i during 
the first Ti units of time because any flight that arrives in the destination center i has to spend a minimum 
of Ti units of time flying in center i before landing. 

Constraint (13) means that there are no incoming flows or departures in a center i during the last Ti 
units of time, because otherwise these flights will still be in air at the end of time horizon (T) in the study, 
due to the minimum flight time requirement. 

In constraint (14), 

N 

^ v fiij (fc)ffj(fc) ~\~ di{k) 

3 = 1 

represents the number of “loss” flights, which equal the number of outgoing flights (X)jli Aj (k) x i(k)) and 
landing flights di(k) in center i at time k. Consequently, the term 

t+Ti / N 

£ ^2/3ij(k)xi(k)+ai(k) 

k=Ti \j=l 

represents the cumulative “loss” flights from time Ti to t + Ti in center i. This is also the cumulative 
“loss” flights from time 0 to t + Ti because during the first Ti units of time there are no loss flights (by 
constraint (12)). As a counter part, the term 

t ( N 

fc=0 \j=l 

represents the cumulative “gain” flights from time 0 to time t. Altogether, the constraint (14) can be 
interpreted as: in center z, the cumulative number of “loss” flights cannot exceed the cumulative number 
of “gain” flights from Ti units of time earlier. The “gain” flights have to spend at least Ti time in center i 
before getting lost. 

In the current study, it is assumed that the minimum flight time is identical for every flight despite the 
fact that flights use different routes that need different flight times, that there are flights taking off and 
landing in the same center, and that flights taking off and landing have different flight times than en route 
flights, etc. However, the proposed mathematical formulation (12)-(14) for minimum flight time can be 
extended to include the above detailed scenarios. For example, using A to denote the starting point of a 
route in center z, and using B to denote the ending point as follows: 


A 


the minimum flight travel time along the route A — > B can be formulated as: 

£ p B j ( k)x AB {k) = 0, fc € {0, • • • , Tab - 1} (17) 

ieiv® 

£ 0 jA (k) X j (k) = 0, k € {T - Tab + 1, • • • , T} (18) 

j €N* 

t+TAB t 

£ £ Mk)x AB (k ) < £ £ PjAikfoik), t € {0, • • • , T - Tab}, (19) 

k=T AB j£N B k = 0 j£N A 
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where XAB(k) is the aircraft count on route A — ► B at time k , is the set of downstream flows connecting 

AT" 4 is the set of upstream flows connecting A. Tab is the minimum flight time on route A — ► 5. 

The minimum flight time is computed using historical air traffic data. Figure 8 shows an example of the 
flight times in the Los Angeles center (ZLA). In this example, the 25-th percentile flight time (20 minutes) 
is chosen as the minimum flight time for all flights through ZLA. Figure 9 shows examples of flight times for 
four different routes inside the ZLA center. 


xIO 


ZLA 



Figure 8. Flight times for all flights in ZLA. Data collected from March, May and September of 2005. 
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Time bin (minutes) Time bin (minutes) 


Figure 9. Upper left: Flight times for flights from ZAB through ZLA to ZLC. Upper right: Flight times for 
flights taking off and landing in ZLA (flights did not exit ZLA). Lower left: Flight times for flights from ZAB 
landing in ZLA. Lower right: Flight times for flights taking off in ZLA and entering ZAB. Data collected from 
March, May and September of 2005. 
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